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Abstract 

The non-planar Feynman diagram with seven massless, scalar propagators and 
four on-shell legs (the crossed double box) is calculated analytically in dimensional 
regularization. The non-planar diagram with six propagators is also discussed. 

1. Introduction 

Knowledge of two-loop massless Feynman diagrams with four massless external legs is 
one of the key ingredients required for a next-to-next-to-leading order calculation of 2 jet 
production rates in hadron-hadron collisions 0, as well as for other processes, including 
Bhabha scattering in the high energy limit. Among those diagrams, the most compli- 
cated are the planar and non-planar double box diagrams, both of which contain seven 
propagators. Recently, significant progress in this field has been made by Smirnov, who 
has derived an analytical formula in terms of polylogarithms for the planar double box 
(see also ref. ||). More simple planar diagrams, with less than seven propagators, are 
discussed in refs. [|], [5]. Results for non-planar diagrams have until now been missing. 





The ultraviolet divergences of the planar and non-planar double box diagrams in space- 
time dimensions d — 7, d — 9 and d — 11 were calculated in ref. 0. In d = 4, the diagrams 
are ultraviolet finite, but they have a complicated structure of collinear and infrared di- 
vergences. In dimensional regularization, with d = 4 — 2e, the divergences appear as poles 
1/eP, with j = 1,2,3,4. As has been demonstrated by Smirnov in the planar case, the 
method of Mellin-Barnes contour integrals is a very natural and convenient technique to 
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isolate these singular contributions. The main purpose of this paper is to calculate the 
non-planar seven propagator diagram, 

b np = f f id. d d, 1 /,\ 

7 J J (k + i + Pl + P2 y(k + 1 + P2 y(k + iy(k - P3 ) 2 k^i - Pi yp ' 1 ) 

using a similar Mellin-Barnes approach. The external momenta pi . . . p± are lightlike 
ipl — 0) an d we use Mandelstam variables defined by s = (pi + p 2 ) 2 , t — (pi + P3) 2 and 
u = (px +Pa) 2 . 

Starting from a Mellin-Barnes representation, which we derive in section 2, we obtain 
analytical formulas in terms of logarithms and polylogarithms for the seven propagator 
diagram B^ p , which we present in section 5, by a series of steps that are explained in 
sections 3 and 4. Then, in section 6, we briefly discuss the non-planar diagram with six 
propagators. Some final comments are made in section 7. 

2. A Mellin-Barnes representation 

Our main tool will be the Mellin-Barnes representation of a power of a sum as a contour 
integral, 

1 r°° T(-a)T(is + a) 

where the integration contour separates the poles of r(— a) from those of T{v + a), and 
v4i t 2 are complex numbers such that | arg(v4i) — arg(A2)| < n. The Mellin-Barnes repre- 
sentation of a sum of several terms, (A± + . . . + A n )~ u , is obtained easily by iteration of 
the basic formula (0). Another formula that we will need is Barnes's first lemma (see, e.g. 



I rioo 

— / da r(oi + a) T(o 2 + a) T(a 3 - a) T(a 4 - a) 

Z7lt J-ioo 

r(o"i + cr 3 ) T(a 1 + 0-4) r(fx 2 + 0-3) r(cr 2 + a 4 ) 



r(ai + a 2 + 0-3 + a 4) 



(3) 



Again, the contour should separate the increasing series of poles (of r(o"3 j 4 — a)) from the 
decreasing ones (of r(o"i j2 + a)). 

Introducing Feynman parameters in the standard way, the non-planar double box 
diagram can be written as 

B" p = (-l)V d / 2 ) 2 r(7-d) r dx 1 ...dx 7 6(l-x 1 ...-x 7 )D 7 -'fc d ' 7 (4) 

Jo 

D = (xi + x 2 + x 3 )(x 4 + x 5 + x 6 + x 7 ) + (x 4 + x 5 )(x 6 + x 7 ) (5) 

C = —tx 2 X$XQ — UX 2 X4X 7 — SX1X5X7 — SX3X4X6 

— S XiX 3 (x 4 + X 5 + Xq + x 7 ) . (6) 

According to the causal ie-prescription, a small imaginary part should be added to C . 
We do this by giving s, t and u positive imaginary parts. At this stage, we consider s, 
t and u to be three independent variables. At the end of the calculation, after we have 
taken the limit where s, t and u become real, one variable can be eliminated by using the 
relationship s + t + u = 0. 
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The next step is to represent the factor C d 7 as a product of five more simple factors, 
at the expense of introducing four Mellin-Barnes parameters: 



T(7-d)C d - 7 = j-— I dM d£ 2 / da I dp 



too ptoo rioo ptoo 



too J — too J —too 



xr(-6) r(-6) r(-a) v{-p) r(7 - d + & + £ 2 + a + p) 

X [—t'X2'Xf > X%[ \-UX2XiX-jf 2 [-S XiX^Xj] 01 [-S X3X4XQ] 

X [— S XiX 3 (x 4 + X 5 + X 6 + x^^-^-^-a-P ^ ^ 



Now the Feynman parameter integrals can be done in terms of T-functions, and we obtain 

N 
7 

ir d T 2 (-e) 



the following Mellin-Barnes representation for B NV ■ 



7 r(-2e) - 3e) 



too rioo 



K = Tw-U I d£xl d& / da I dp 



—too J —too J —too 



(2m) 

x (- S )- 3 - 2 ^ 2 (-tf (-uf r i r 2r3r 4 r 5 r 6 rg 8 r 9 r 10 r 11 r 12 r 1 3 ^ (g) 
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where 



ri = T(3 + 2e + a + P + & + 6) T 8 = r(-2 - 2e - a - 6 - 6) 

r 2 = r(2 + e + a + p + Ci + 6) r 9 = r(-2 - 2e - /? - a - 6) 

r 3 = r(i + £i + 6) r 10 = r(-a) 

r 4 = r(i + a + ^) r n = r(-/5) (10) 

r 5 = r(i + a + 6) r 12 = r(-6) 

r 6 = r(i + /5 + ei) r 13 = r(-e 2 ) 

r 7 = r(i + /3 + £ 2 ) r 14 = r(2 + a + /? + ei + 6). 

In order to keep things as simple as possible, we choose the contour for each of the 
variables £1, £ 2 , a and /? to be a straight line parallel to the imaginary axis, with a real 
part that does not depend on any of the other variables. Eq. (|9|) is certainly valid if the 
real parts of the contours are chosen in such a way that the real parts of the arguments 
of r\ . . . r 13 are all positive. It is possible to satisfy these conditions simultaneously when 
-3/2 < 9fc(e) < -1/3. For example, if e « -1/2, choosing = = -5/12, 

dt(a) = dt(P) = —1/4, satifies all the conditions. Obviously, there are infinitely many 
other, equally suitable possibilities, but for the sake of definiteness, we will stick to these 
contours. 

3. Analytic continuation in e 

Having derived the Mellin-Barnes representation (g) near e = — 1/2, we must now perform 
an analytic continuation in e to the region around e = 0. To do so, we follow the poles of 
the T-functions as we increase e. Whenever a pole crosses one of the integration contours, 
which we keep fixed, we add its residue in the corresponding variable to the right hand 
side of eq. @. This residue term still contains poles in its remaining variables which can 
cross other integration contours as e is increased further. They are treated in the same 
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way. We end up with a large collection of residue terms, plus the original fourfold integral. 
Some of the residue terms contain factors that diverge at e = 0. However, these singular 
factors do not depend on £1, £ 2 , ot or (3. The remaining integrals can simply be expanded 
in e under the integral sign because our contours do not go through any poles at e = 0. 
(There are never problems with convergence at the limits £i,£ 2 ,a,/3 — > ±ioo). 

Below, we will give some examples of the residue terms that appear. But first, let us 
see what happens when we substitute e = in the integral (|9]). The variables a and (3 
are coupled via r 1; T 2 and T 14 , but these T-functions collapse into 

r 1 r 2 /r? 4 = 2 + a + /3 + & + 6 = (1 + « + &) + (1+0 + 6). (11) 

Consider the first term. The factor (1 + a + £1) can be absorbed into T 4 , and then both 
the a and the f3 integrals can be performed using Barnes's first lemma. Moreover, the f3 
integral gives zero: 

I MOO 

— / d/3r(i + /3 + e 1 )r(i + /3 + e 2 )r(-2-/3-ei-6)r(-/3) 

r(-i - &) r(-i - 6) r(i + e0 r(i + 6) n ri9 . 

r(o) " u - 1 j 

(Here, we are taking a contour that separates the increasing and decreasing series of poles 
as required by Barnes's lemma, not a straight line). Analogously, the second term on the 
right hand side of eq. ([TT|) vanishes when integrated over a. 

This cancellation is very nice because it means that if we only need the non-planar 
box up to the constant term in e, once we have done the a and (3 integrals, we do not 
need to calculate any double integrals of the form 

d£i r d£ 2 (s)- 3 -^ (-t)^ (~uf .... (13) 

J — ioo 

Only terms in which a residue has been taken in either £i or £ 2 , or in both, are ever needed. 
Such terms, whose dependence on at least one of the variables s, t and u is trivial, are 
much more simple than the double integrals (|T^). 

We now return to the task of collecting the various residue contributions. The first 
poles to cross contours are at a = — 2 — 2e — £i — £ 2 and (3 = — 2 — 2e — £i — £ 2 , coming from 
T 8 and T 9 , respectively. (The poles of r x and T 2 are harmless because they depend on e 
with the opposite sign, and therefore move away from the contours when dt(e) increases). 
To account for the residues of these two poles, we replace K — > K + K$ + Kg + K 8 g, where 

(14) 



(27T?) 3 J-ioo 



X 



r 2 

1 14 



=-2-2e~5i-6 

Kg is a similar term, where the residue in (3 is taken, and 

I moo rioo 

(2^2)^ J-ioo 



Tir2r3r4r5r6r7 riornri 2 ri3 

X 



r 2 

1 14 



(15) 



«=/?=-2-2e-a-6 
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Let us concentrate on if 8 . Note that in eliminating a, we have changed the e-dependence 
of the remaining T-functions. As a consequence, the first poles of T 2 = T(j3 — e), r 4 = 
T(— 1 — 2e — £ 2 ) an d T 5 = T(— 1 — 2e — £1) cross the /?, £2 and £1 contours, respectively. 
These crossings are taken into account by replacing if 8 — > K 8 + K 82 + K M + if 85 + 
K 8i2 + -^852 + -^845 + -^8452- In a number of these terms, there are still more poles that 
cross contours before e finally reaches 0. They are dealt with by making the following 
further replacements: if 82 -> if 82 + if 829 ', if 84 -> if 84 + if 843 , if 842 -> if 842 + if 8423 , 
if 85 -> if 85 + if 853, if 8 52 -> if 852 + ^8523, and if 845 -> if 845 - if 8459 - In each term, the 
subscripts indicate the T-functions of which the residue of the first pole has been taken. 
An exception to this is if 829 /, where the prime means that we have taken the residue of 
the second pole of T 9 (i.e., where the argument equals -1) instead of the first pole. The 
last term, if 8459, gets a relative minus sign because the pole in T 9 crosses the (3 contour 
backwards.0 

Due to the symmetry of the problem in a and (3, there is a one to one correspondence 
between the descendants of if 9 and those of if 8 , and their contributions to the final result 
are the same. The analysis of if 8 g runs along similar lines. Here, some poles are second 
order, with residues involving derivatives of T-functions. 

4. Evaluating the if 's 

In all, we find 39 contributions. (This number depends on the contours one uses. Had 
we chosen different ones, there would have been other poles crossing contours and we 
would have obtained a different, but equivalent, collection of residue terms). The simplest 
contributions are those where a residue in four variables has been taken, e.g. 



if 845 2 = (sy 3 - 2 ^-^ (-tf (-uf ™ WloWl3 



r 2 

1 14 



£1 = & = -1 - 2e 

a = 2e, (3 — e 



= {-sy 1+2e {-ty 1 - 26 i-uy 1 -' 6 r(i + e ) r(-i - 4e) r(e) r(-2e) r(-e) r 2 (i + 2e) . (16) 

Note that in this particular case, the e expansion starts off with a 1/e 4 pole. 

In several terms, there is one integral left. For instance, if 84 3, which we get by taking 
residues at a — — 2 — 2e — £1 — £ 2 > £2 = — 1 _ 2e and £1 = 2e, is given by 

^843 = (sy 2 - 2e (-t) 2e (-uy 1 - 2 " r(-l - 4e) T 2 (l + 2e) T(-2e) 
x J_ r d3 r(l + (3) T((3 - e) T(l + 2e + $) T(-l - 2e - 0) T(-/3) 
2mJ-ioo P r(/3-2e) ' 1 ' 



To calculate the (3 integral, we expand the integrand up to the second order in e. Then, 
closing the contour to the left or right and summing the residues of the poles inside it, 
we obtain harmonic series which can be expressed in terms of C(2), C(3) and C(4)- 

An example in which the final integral still involves powers of Mandelstam variables 

is 



27riuJ-ioo £i(l + £ 



where we have already set e = 0. This integral can also be calculated by closing the 
contour, and the result can be expressed in terms of polylogarithms Li n M and Nielsen's 



Eventually, i^8459 cancels against a descendant of K t 



89 ■ 
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generalized polylogarithms S n>p |9], |10| , which are defined by 



s^jz^/' ^'M^t'-*) . (19) 

(n — jo t 

The cases Sx,2, Si^ and 82,2 which occur in this paper can all be written as combinations 
of Li 2 , Li 3 and Li4. 

In the term K%±, we perform the first integration by Barnes's lemma (in a slightly 
modified form, see [Q]): 

Km = 1 r d6 (_ s) — ( _ t)6 r»(i + &)!-(-&) 



(2tu) 2 u 7-ioo " v ' v ' 6(1 + 6) 

/ioo 
d(3T(i + p)r(i + p + 6) r(-i - p - a) r(-/3) 
-ioo 

1 I*" d6 (-.)-* Efa + f4g(z6) w + 6) - m) . (20) 



2m u J -ioo + 

The £\ integral can be then be done by closing the contour, as in the previous example. 
Three other K's are related to K M by a <-> (3 and/or £x <-> £ 2 . 

Finally, there are four termsQ in which both the £1 and the £ 2 integrals still survive, 
namely K, Kg, Kg and K$g. However, by the mechanism explained above (see eq. (|12D), 
these terms cancel against each other at e = after integration over a and (3. 

5. Result for the seven propagator diagram 

Adding together all the contributions discussed in the previous section, we obtain an 
expression for B^ p , up to the constant term in e, in terms of (poly) logarithms of (ratios 
of) s, t and u, which are, at this point, still independent variables. In particular, the 
result is also valid when s, t and u are all negative. This region is, of course, unphysical, 
but it is nevertheless interesting to consider because it provides us with an opportunity 
to check our result against a calculation by Binoth and Heinrich jD]]. They have obtained 
results for the negative s, t, u region by extracting the divergences from the Feynman 
parametric integral (H) and then performing a multi-dimensional numerical integration. 
At the symmetric point, s = t = u = — 1, we find 

, w „ 1, J 7 3 1 (-1 47 2 \ 1 /105 37 2 89 

?i VP = 7T r Z 4-9f 1 7T 2 14- 1 7T 2 1 



- ' ar < 3 + 2 ')(4?-? + ?l--^J + 7l- + T^-T^ 

2 6V 7 2 3 288 4 sw 

+ y vr 2 log(2) - y vr 2 log 2 (2) + H i og 4 (2 ) + 22 Li 4 (1/2)} 

, , , f 1.75 3 22.828 113.63 1 , , 

= 7r d r(3 + 2e)|-^-- — + — 395.26| . (21) 

The coefficients in ( ^T|) agree with the results obtained by the numerical method [111]] . We 
also find good agreement at the asymmetric point (s,t,u) = (—1, —2, —3). In both cases 

3 We are leaving out K S2 and K 92 , which contain a factor of T(— e) in their denominator and are 
therefore already of order e. 
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the precision of the numerical method ranges from about 1 per mille for the l/e 4 -term to 
1 per cent for the constant term. 

Since the complete formula for general s, t and u is very long and not relevant for 
physical applications anyway, we do not present it here. Instead, we specialize to the 
physically relevant case where s + 1 + u = 0. We also use transformation formulas for the 
arguments of the polylogarithms in order to avoid ending up on a cut when we make s, 
t, and u real. Depending on which legs correspond to incoming and outgoing particles, 
there are three different kinematical regions to consider: (i), where u, t < < s; (ii), 
where u, s < < t; and (iii), where s, t < < u. In all regions, we have 

B ~P = (^)V (1+£) {jL + J|}. (22) 

We will use the abbreviations T = log(— t/s), U = log(— u/s), V = log(— u/t). 
In region (i), with u, t < 0, and s = —t — u, F t is given by 

F t = S - 2e {^ + ^(5T + 7f/) + i(6f/-T 2 -4Tf/-f/ 2 -^7r 2 ) 

+ - (-24 U - 12 TU - ^T 3 + 3 TU 2 - U 3 + ^tt 2 (5 T - 29 U) 

19 

-2 T Li 2 (-t/s) + 2 Li 3 (-t/s) - 2 S lj2 (-t/s) + —((3) 
+96 U + 48 TU - 12 T 2 U + 12 TU 2 - 4 U 3 + ^T 4 + ^T 3 U - T 2 U 2 

-(13T + 33C/) C(3) + (-48T+ 17T 2 + 12TC/+ IOtt 2 ) U 2 (-t/s) 
+ (48 - 60 T - 12 17) Li 3 (-t/s) + (-24 + 28 T - 6 [/) Si, 2 (-t/s) 
+86 Li 4 (-t/s) - 26 Si, 3 (~t/s) - 36 S 2 , 2 (-t/s) 



2 1 1 / 31 

- + — (6 - T + f/) - - 24+12 T + 2T 2 + 2 Tf/ + 2 rj 2 + — tt : 

e d e V 6 

+96 + 48T+ 12T 2 -24Tf/+ — T 3 + 2U 3 

3 

+7T 2 (-2 + yT - + (-24 + 14 T + 18 U) Li 2 (-t/s) 
-32 Li 3 (-t/s) + 44 Si, 2 (-t/s) - 89 C(3) 



(23) 



and F u can be obtained by interchanging (t <-> -u) in the expression for F t . 
In region (ii), with tt, s < 0, and t = —s — u, F t and F u are given by 

F t = t~ 2e {^ + ^(2T+^V^+^(6T + 6V + 2T 2 + TV-V 2 + ^7r 2 

+ i (-24 T — 24 V — ^T 3 — 2 T 2 V - 2 TV 2 - V 3 - ^tt 2 (23 T + 41V) 

15 



+2S li2 (-s/t) + y C(3)) 
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+96 T + 96 V - 4 T 3 - 12 T 2 V - 4 V 3 - — T 4 - — TV + T 2 \/ 2 + 2 TV 3 + -V 4 

6 3 3 

+7T 2 (-14 T - 14 V - — T 2 - —TV + -V 2 ) - — 7T 4 
V 3 3 6 / 120 

- (24 + 45 T + 39 V) C(3) + (24 T + 6 T 2 - 18 TV + 13tt 2 ) Li 2 (-s/t) 

+ (24 + 12 T - 18 V) Li 3 (-s/t) + (24 - 44 T + 6 V) Si )2 (-s/t) 

+ 12 Li 4 (-s/t) + 26 Si, 3 (-s/t) - 62 S 2 , 2 (s/t) 



54< r+4V, » 



+i (l2 T + 12 V + 4 T 2 + 2 TV - 3 V 2 + ^tt 2 - 2 Li 2 (-s/t)) 

-48T-48V + 12T 2 + 24TV- 12 V 2 + ^T 3 + 2 TV + -V 3 
+7T 2 (4 - 6V) + (-48 + 14 T + 12 V) Li 2 (-s/t) 



+32 Li 3 (-s/t) + 28 Si, 2 (-s/t) - 15 C(3) 



} , (24) 



—2 1 / 5 \ 1 / 31 
— + -(2T+-V)+ 72 (QT + 2T 2 -TV-V 2 + -', 



+- (-24 T - 12 TV - -T 3 - 4 TV - 2 TV 2 - ^ V 3 - -n 2 (23 T + 31 V) 
6 V 3 3 6 

+2 S 1>2 (-*/*) + yC(3)) 

11 5 10 2 

+96 T + 48 TV - 4 T 3 + 12 TV 2 T 4 - -TV + 4 T 2 V 2 + — TV 3 + -V 4 

6 3 3 3 

/ 92 7 31 \ 311 

+7T 2 -14 T - — T 2 - -TV + —V 2 - — 7T 4 - (24 + 45 T + 21 V) C(3) 
V 3 3 by 120 

+ (-24T- 6T 2 - 18TV + llvr 2 ) Li 2 (-s/t) - (24 + 12T+ 18 V) Li 3 (-s/t) 
+ (48 - 32 T + 26 V) Si, 2 (-s/t) - 12 Li 4 (-s/t) + 86 Si, 3 (-s/t) - 50 S 2j2 (-s/t) 

■Htt g + I(- 6 -T + 4V) 

+- ^24 + 12 V + 2 T 2 + 2 TV - V 2 + — 7T 2 - 2 Li 2 (-s/t) 

-96 - 48 V - 12 V 2 + -T 3 - 4 TV - 2 TV 2 - - V 3 

3 3 

+7T 2 (6 + 6 T - 8 V) + (-24 + 26 T - 8 V) Li 2 (-s/t) 



+44Li 3 (-s/t) - 24S 1)2 (-s/t) + 75C(3) 



(25) 



The formulas for F t and F u in region (iii) can be obtained easily from the ones for region 
(ii) by using the fact the B^ p is symmetric under (t <-> u). This follows from the symmetry 
between the legs of the non-planar double box carrying momenta p 3 and p&. 

6. Six propagator diagram 



8 



The non-planar diagram with six propagators, 



B 



NP 



d d k d d l 



(k + l + p 2 ) 2 (k + lY{k - p 3 ) 2 k 2 (l - Va) 2 1 2 ' 



(26) 



is far less difficult to calculate than B^ p . Perhaps the easiest way is by considering a 
scale transformation, 



B^ p (X Pl ,Xp 2 ,Xp 3 ,Xp,) = A 2(d_6) Bq P (pi,P2,P3,Pi) 



(27) 



Differentiating both sides of eq. (^7]) with respect to A at A = 1 gives an identity which 
can be used to reduce Bq P to a sum of three planar diagrams: 





(28) 

In the diagrams on the right hand side of eq. fl2"8|), the propagator on the diagonal is 
squared. Such boxes with a diagonal propagator are calculated in ref. H. They are, in 



fact, rather similar to one-loop massless box diagrams |L2|. The pole at e = —1/4 is 
a reflection of the linear infrared divergence?] in (^) coming from the region where the 
loop momenta k and / are both soft. Another curious feature of eq. (^) is that although 
the three terms on the right hand side contain third order poles at e = 0, the leading 
singularity of the sum is only 1/e 2 . 

The six propagator diagram is completely symmetric under permutations of its exter- 
nal momenta. Therefore, it is sufficient to consider only region (i), where u, t < < s and 
s = —t — u. We obtain 



r>NP 
B 6 



d/2 



2 T(1 + 2e)r 3 (l -e) 3 ( s~ 



(T + iff) (1 — 2eU) 

r(l — 3e)(l + 4e) t l [ su V 
+2e 2 (^2 Li 3 (-*/*) - 4 Si,2 (-*/*) + 2 C(3) - ^T 3 

2U 2 {-t/s) + U 2 



+T [-2 Li 2 (-t/s) + U 2 - 7T 2 



VK 



-TT 



(29) 



The reduction formula (f2H) can be generalized to the case where the six propagators 



in eq. ( |2q ) are raised to arbitrary integer powers. 
7. Discussion 

The analytical property that distinguishes the non-planar diagrams from the planar ones 
considered in refs. [^|, |||, [| is the fact that they have cuts in three channels, s, t, and u, 
rather than just in two. As a consequence, they have imaginary parts in all three physical 
regions, (i), (ii) and (iii). This is why, when we derived the Mellin-Barnes representa- 
tion (S) for the seven propagator diagram, we were led to considering it as a function of 



It would diverge as l/m reg if we used a mass m reg to regularize it. 
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three independent variables. A priori, one would think that going from a two-scale prob- 
lem to a three-scale problem could make it vastly more difficult to solve, but it did not, 



because, thanks to the vanishing of (12), all genuinely three-scale contributions cancelled 
out. 

Generalizing the Mellin-Barnes representation for the seven propagator diagram to 
arbitrary powers of propagators is completely straightforward, and merely amounts to 
shifting the arguments of the T-functions ( |T0| ) by some constants. It is likely that can- 
cellations similar to (|12|) will take place for other integer powers, so that the method 
described in this paper could be used for those well. 

I am very grateful to T. Binoth and G. Heinrich for their help with the numerical 
checks, and for many discussions. I would also like to thank J.J. van der Bij, A.I. Davy- 
dychev, T. Gehrmann, K. Melnikov, E. Remiddi and V.A. Smirnov for helpful discussions 
and correspondence. This work was supported by the DFG-Forschergruppe "Quanten- 
feldtheorie, Computeralgebra und Monte-Carlo-Simulation" . 
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